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Abstract 

Four results are given that address the existence, ambiguities and con- 
struction of a classical R-matrix given a Lax pair. They enable the uniform 
construction of R- matrices in terms of any generalized inverse of adL. For 
generic L a generalized inverse (and indeed the Moore- Penrose inverse) is 
explicitly constructed. The R-matrices are in general momentum depen- 
dent and dynamical. The construction applies equally to Lax matrices 
with spectral parameter. 



1 Introduction 



The modern approach to completely integrable systems is in terms of Lax pairs 
L, M and R-matrices. Here the consistency of the matrix equation L = [L, M] 
expresses the equations of motion of the system under consideration. The great 
merit of this approach is that it provides a unified framework for treating the 
many disparate completely integrable systems known. Given a 2n-dimensional 



phase space Liouville's theorem |L2], [LJ, which ensures the existence of action- 
angle variables, requires that we have n independent conserved quantities in 
involution; that is they mutually Poisson commute. As a consequence of the 
Lax equation the traces Tr L k are conserved and these are natural candidates for 
the action variables of Liouville's theorem. (In practice the action variables are 
typically transcendental functions of these traces.) It remains however to verify 
that these traces provide enough independent quantities in involution. Verifying 
the number of independent quantities is usually straightforward and the remain- 
ing step is then to show they mutually Poisson commute. The final ingredient of 
the modern approach, the R- matrix |jl5| , guarantees their involution. If L is in a 
representation E of a Lie algebra g (here taken to be semi-simple), the classical 
R-matrix is a E®E valued matrix such that 

[R,L® 1] - [R T ,1®L] = {L®L}. (1) 

(The notation is amplified below.) Then 

{Tie L\ Ti e L m } = Ti E ® E {L k ? L m } = km Ti e ®e L k ~ l ® L m ~ l {L ® L} = 0, 

which vanishes due to the cyclicity of the trace. By a result of Babelon and 
Viallet H|, such an R-matrix is guaranteed to exist if the eigenvalues of L are 
in involution. The Liouville integrability of a system represented by a Lax pair 
has been reduced then to finding any solution to ([!]) and counting the number 
of independent traces. Further, the R-matrix is an essential ingredient when 



examining the separation of variables of such integrable systems [16|, [TTJ . 

Unfortunately the construction of R-matrices has hitherto been somewhat 
of an arcane art and many have been obtained in a case by case manner ||. 
The purpose of this note is to present four results that address the existence 
and construction of solutions of (|1|) and hence the Liouville integrability of the 
system under consideration. They yield a uniform construction of R-matrices. In 
fact the R-matrices satisfying (JI|) are by no means unique and our construction 
characterises this ambiguity. The approach applies equally to R-matrices with 
spectral parameter. We will illustrate these results with a simple example. At the 
outset we remark that the R-matrix solutions to (|]) are generically momentum 
dependent. Within this family of solutions some may be particularly simple: 
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they may for example be constant (as in the Toda system ||) or momentum 
independent. We are content here with providing the construction of an R- matrix 
given a Lax matrix L and so answering the question of Liouville integrability: 
we do not seek to further specify the momentum or position dependence of 
the solution. For the elliptic Calogero-Moser model there are in fact || no R- 
matrices that are independent of both momentum and spectral parameter (for 
more than 4 particles) and this illustrates the fact that simple assumptions on 
the parameter dependence of an R-matrix need not be natural. Elsewhere we 
will apply these results to the elliptic Calogero-Moser models without spectral 
parameter. 

Our approach is as follows. First we rewrite (|IJ) in the form of the matrix 
equation 

A T X - X T A = B. (2) 

Here A is built out of L and the Lie algebra, the unknown matrix X being solved 
for is essentially the R-matrix in a given basis and B represents the right-hand 
side of (p. Our first result is to give necessary and sufficient conditions for (0) to 
admit solutions together with its general solution. This general solution encodes 
the possible ambiguities of the R-matrix. Because A is (in general) singular our 
solution is in terms of a generalized inverse G satisfying 

AG A = A and GAG = G. (3) 

Such a generalized inverse always exists. (Accounts of generalized inverses may 
be found in || [13], @-) Indeed the Moore- Penrose inverse -which is unique 
and always exists- further satisfies (AG)* 1 = AG, {GAY = GA. Observe that 
given a G satisfying (H) we have at hand projection operators P\ = GA and 
P2 = AG which satisfy 

APi = P 2 A = A, P 1 G = GP 2 = G. (4) 

Our second result shows that the choice of generalized inverse G only alters the 
R-matrix within the ambiguities specified by the general solution, and so any 
generalized inverse suffices to solve (Q) and hence construct an R-matrix. At 
this stage we have reduced the problem of constructing an R-matrix to that of 
constructing a generalized inverse G and our third result constructs such for a 
generic element L of g. Because the Moore-Penrose inverse is unique, our fourth 
result is to present this inverse for generic L though we shall not need to use this 
in our application. 

The letter is organised as follows. In the next section we present the four 
results. The proofs of the first two are somewhat lengthy and algebraic and will 
be presented elsewhere [[F); the proofs of the remaining two are easier to outline. 
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In section 3 we apply these to give the R-matrix for generic L. In section 4 we 
extend the results to include a spectral parameter. Section 5 is an illustrative 
example. We conclude with a brief discussion. 

2 Four Results 

Our first task is to identify ([[]) with (|2]). Let T M denote a basis for the (finite 
dimensional) Lie algebra g with [T M , T v ] = T A defining the structure constants 
of 0. Set </>(T M ) = X^, where <p yields the representation E of the Lie algebra 
0; we may take this to be a faithful representation. With L = L^X^ the 
left-hand side of ([I]) becomes 

{L® L} = Y,{V t ,L"}X (t ®X v 

while upon setting R = R^X^ ® X u and R T = R U ^X^ <g> X v the right-hand side 
yields 

[R, L ® 1] — [R T , 1 ® L] = R^{[X ll ,L\®X l/ -X v ®[X lM ,L\) 

= R^L\[X^X X }®X U -X U ®[X^X X ]) 
= (R^c^-R^L^X^X,. 

By identifying A"" = c u ^ x L x = -ad(L) u ^ B» v = L u } and X^ = R^ u we see 
that (1]) is an example of (0). 

Having shown how to identify ([!]) with the matrix equation @ we may now 
state our first result. 

Result 1 The matrix equation (^) has solutions if and only if 

(CI) B T = -B, 

(C2) (1 - P?)B(1 - P x ) = 0, 

in which case the general solution is 

X = \G T BP 1 + G T B(1 - Px) + (1 - Pl)Y + (P?ZP 2 )A (5) 

where Y is arbitrary and Z is only constrained by the requirement that P 2 T ZP 2 
be symmetric. 

Although the general solution appears to depend on the generalized inverse G 
we in fact find 
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Result 2 If G is any other solution of (||) with attendant projection operators 
Pi 2 then (fjj may also be written 



1 
2 

where 



X = ^G T BP 1 + G T B(1 - Pi) + (1 - P 2 T )F + P?ZP 2 A 



Y = (1 - P 2 T )F + P[ZP 2 A + G T P(1 - -Pi) Z = Z + ^(G T BG - G T BG). 

Thus Z is again symmetric and we have a solution of the form (|5]). 

In the R-matrix context the matrix B is manifestly antisymmetric because 
of the antisymmetry of the Poisson bracket and so (CI) is clearly satisfied. We 
have thus reduced the existence of an R-matrix to the single consistency equation 
(C2) and the construction of a generalized inverse to ad(L). We turn now to the 
construction of the generalized inverse. 

Let X^ denote a Cartan- Weyl basis for the Lie algebra jj. That is {X^} = 
{Hi, E a }, where {Hi} is a basis for the Cartan subalgebra f) and {E a } is the set 
of step operators (labelled by the root system $ of q). The structure constants 
are found from 



[Hi, E a ] = ai E a , [E a , E. a ] = a y -H and [E a , E p ] = N afi E a+p if a+(3 G $. 

a+fi 
af3 



Here N a a = c*~\P . With these definitions we then have that 
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(6) 



(adL) 



a V -ajL a A 



1 




^ 


u T \ 




= ( 


v v 


A ) 



where we index the rows and columns first by the Cartan subalgebra basis 
{i,j : l...rank$j} then the root system {a, (3 G $}. We will use this block 
decomposition of matrices throughout. Here u and v are |$| x rankg matrices 
and we have introduced the |$| x |$| matrix 

A a p = a-L5% + c a a _^L a -e, (7) 
where a ■ L — Yli=i° a i-L l - With these definitions we have 

Result 3 For generic L the matrix A is invertible and a generalised inverse of 
adL is given by 

1 \f \ / 1 -u T k- 1 \ _ f 
-A" 1 *; 1 j I A 1 Mo 1 j ~ 1 A" 1 
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We establish the result by first showing that for generic L 

<-) - ( I T ) ( Si ) ( A, : ) = m 

it then follows that (§) is a generalised inverse for adL by direct multiplication. 
Now for any matrices m and A we have the general factorisation fH| |5|] 

171 uT \ - ( 1 u T E\('m-u T Ev u T {l - HA) \ / 1 
v A ) ~ V 1 J V (1 - AE)v A J \ H v 1 

where 5 is a generalised inverse of A. In particular, when m = and A is 
invertible (and so H = A -1 ) this shows that 

u T A / 1 u T A- x "\ / —u T hr l v \ / 1 0\ ^ 



v A J \0 1 J \ A J \ A^v 1 

Thus (0) and hence the result follow by establishing that A is generically invert- 
ible and that 

u T A~ l v = 0. (11) 

From (^) we see that A is the perturbation of a diagonal matrix and so is 
generically invertible: the zero locus det A = is a polynomial in the coefficients 
of adL and so the complement of this set is dense and open. For such an invertible 
A we thus have 

rank A = dim A = dim g — rank g. (12) 
Now the maximum rankQ of the matrix adL is dim g — rank g [I7| . From ( ID ) 



we see that rank A + rank(w T A -1 i;) = rank(adL) < dimjj — rankjj and so from 
fli~2l) we deduce that rank A = rank(arfL). Therefore rank(u T A -1 v) = and 
consequently ( |TTJ ) must hold. The result then follows. 

An alternate factorisation of adL is possible for the generic L under consid- 
eration. Utilising ( |TTD we find that 



(adL) 



'^(A^v 1)=EAF. 



Employing a result of MacDuffee || this full rank factorisation then yields: 
Result 4 For generic L the Moore-Penrose inverse of adL is given by 

F^(F F j! )~ 1 A~ 1 (E j ' E)~ 1 E^ 

( u T A~ x \ 

where E = ( J and F = ( A^v 1 ) . 

1 If det(t — adL) = X]j™ B Pj (L) & is the characteristic polynomial of adL, the regular semi- 
simple elements of a semi-simple Lie algebra g are those elements for which p lan ^ g (L) ^ 0. 
These elements are also of rank dimg — rankg and form an open dense set in g, but this 
condition is different from det A ^ 0. 
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3 The R-matrix 

We now bring together the results of the previous section to present the R- 
matrix for a generic L when this exists. From the fact that A = —(adL) T 
a generalized inverse of A is given by minus the transpose of the generalized 
inverse (H). Utilising our earlier notation this means that we have the projectors 



, \ / v T A~^ 

1 1 ( A- 1T u 1 ) ' P * ~ \ 1 

Let us express the Poisson brackets of the entries of L in the same block form in 
the Cartan-Weyl basis: 



B= f C f ) -D 



T 



where B a ^ = {L a ,L J } = fi a j and so on. The constraint (CI) is manifestly 
satisfied. 

The constraint (C2) is now (the rankg x rankg matrix equation) 

((72) = (+/A- 1t m-/A-V+/A"VA^ 1t «. (13) 

Each term in this equation is known and so the equality may be readily checked. 

Supposing the constraint ((72) is satisfied we then find from (|5|) the general 
R-matrix takes the form 

R=( ° ° W P q 

\ -A-V + !A~VA -1T w -|A~V J V -A - V - Fu —A~ l vq — FA T 

(14) 

The second term characterises the ambiguity in R where we have parameterised 

j and Z = \ <i ) ' ^ ere ^ ne ma t r ices 
p, g are arbitrary while the entries of Z are such that 

F = A _1 f at> T A~ 1T + d + A- 1 ^ + ct; T A- 1T (15) 

is symmetric. 



4 Inclusion of Spectral Parameter 

For simplicity we have presented our construction for Lax pairs with no spec- 
tral parameter but it is straightforward to incorporate such a parameter. The 
relevant equation to be solved for is now 

{L{u) ® L{v)} = [R(u, v), L{u) (8) 1] - [RT(u, v), 1 ® L(v)}, (16) 
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where if R(u,v) = W JV (u,v)X ll <8> X v then R?(u,v) is defined byf] R?(u,v) = 
R u ^(v,u)X fM (g)X u . Now 

B^(u,v) = {L»(u),L v (v)} = -B^{v,u) 

and because L(u) depends on u alone the generalized inverse now also depends 
on the spectral parameter as G = G(u). The equation we now wish to solve is 

A T {u)X{u, v) - X n (u } v)A(v) = A T (u)X(u, v) - X T (v } u)A(v) = B(u, v), 

and this has the analogous solution 

X(u,v) = \G T {u)B{u,v)P l {v) + G T {u)B(u,v){l- P x {v)) 
+ (1 - P?(u))Y(u,v) + (P 2 T (u)Z(u,v)P 2 (v))A(v) 

if and only if 

(CI') B n = -B, 

(C2') (1 - P[(u))B(u,v)(l - P x {v)) = 0. 

Here Y(u,v) is arbitrary while the symmetry condition now becomes 

(Pl(u)Z(u,v)P 2 (v)r = P 2 T (u)Z(u,v)P 2 (v). 

As in the spectral parameter independent case, this reduces to the requirement 
that 

F(u, v) = A~ 1 (-u)'y('u)a('U, v )v T {v )A~ 1T (t> ) + d(u, v) + A~ 1 (u)v(u)b(u, v) 
+c(u 1 v)v T (v)Ar 1T (v) 

(17) 

be such that F 7 * = F. 

5 An Example 

We conclude with the simple but illustrative example of the harmonic oscillator 
presented as the Lax pair (with spectral parameter) 

^> =((*>("+! ^/t 1 ). "<«>=(?;)■ M 

The consistency of the Lax equation L{u) = [L(u), M(u)] follows from the equa- 
tions of motion of the Hamiltonian H = (p 2 + x 2 + u)/2 = —(u/2) det L(u). 

2 We use 7r to denote both matrix transposition together with the interchange of u and v 
while T denotes ordinary matrix transposition. 
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Although we could equally work with the simple algebra su(2) in this example 
we will take the algebra to be gl(2). Now for any I G gl{2), 



a b 
c d 



aH x + bE l2 + cE 2 i + dH 2 



we find that in our Cartan-Weyl basis 



ad I 



/ 








— c 


b\ 










c 


-b 




-b 


b 


a — d 





V 


c 


— c 





d — a J 



Then 



u 



-c 
b 



-b 

c 



A 



a — d 






d — a 



and I is generic^ provided a — d ^ 0. We note that in the Cartan-Weyl basis the 
permutation operator P = Hi® Hi + Yl a ^ <2) E_ a (which is such that 
P(X <g) F)P = y ® X) takes the form 



/ 1 \ 

10 

1 

\ 1 J 



(19) 



Now for the case at hand b = (p 2 /u) + 1, c = {x 2 /u) + 1, a = —d = ipx/u 
and L is generic for px ^ 0. We may therefore use the expressions computed in 
sections 3 and 4. We calculate (using {p, x} = 1) that 



B(u,v) 





( 








-2ip 2 


2ix 2 


\ 


1 










2ip 2 


-2ix 2 




uv 




2ip 2 


-2ip 2 





Axp 






\ 


-2ix 2 


2ix 2 


— 4xp 





/ 



-fi T {u,v) 
6(u,v) 



and straightforwardly verify that condition (C2') is satisfied. The R-matrix is 
then given by 



R(u, v) 





-A _1 (m)//(w, v) + |A _1 (m)0(m, v) A~ 1T (u)u(v) 



.1A-i 





[u)cj)(u,v) 



+ 



p(u,v) 

— A _1 (m)u(m)p(m, u) — F(u,v)u(v) 
It is regular semi-simple provided (a — rf) 2 + be ^ 



g(tt,u) 

A~ 1 (-u)t>(-u)g('u, u) — -F(m, f )A T (v) 
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The second term again characterises the ambiguity in R and we have param- 
eterised the matrices Y(u,v) and Z(u,v) in an analogous way to the spectral 
parameter independent case of section 3. Substitution of the various quantities 
gives for the first term 



R(u, v) 



/ 











\ 


















— p 2 +v 


— p 2 +v 




i 


i 


\ 


2pxv 
—x 2 +v 


2pxv 
—x 2 +v 


V 


2pxv 


2pxv 


V 



This R-matrix is clearly dynamical. Making use of the the block structure of the 
R-matrix we see that by choosing 



p(u,v) 



-2i 



u 



1 
1 



q(u, v) = 0, F(u, v) 



we arrive at the nondynamical 



R(u, v) 



-2i 



P 



u — v 



U + V 1 



u 



v 2px 



1 
-1 



where P is given by ([191) . 



6 Discussion 

We have presented a uniform construction for a classical R-matrix given a Lax 
pair, thus answering the question of the Liouville integrability of the system in 
terms of the invariants of the matrix L. The method not only gives necessary 
and sufficient conditions for the R-matrix to exist and describes its ambiguities, 
but is algorithmic as well. Given L, first construct adL. Next construct any 
generalized inverse to adL and verify (C2); this is the necessary and sufficient 
condition for an R-matrix to exist: it is given explicitly by (||). Further we have 
given a generalized inverse for generic L in (^); genericity is easily checked by 
evaluating det A ^ 0, where A is the restriction of adL to the root space (given 
by (^)). The ambiguities in the R-matrix have been specified. We remark, that 
the block nature of the R-matrix allows us to easily verify putative ansatz for a 
given R-matrix. 

Thus far our discussion has been limited to linear R-matrices and we briefly 
discuss the application to quadratic r-matrices, ie. the solutions to 

{L® L}= [r,L®L] = [r A ,L®L], (20) 
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where ta = (r — r T )/2. (It follows from the anti-symmetry of {, } that [r + 
r T , L (g) L] = 0.) As discussed in [f|], the quadratic R-matrix calculation may be 
reduced to the linear R-matrix situation. In particular the R-matrix 

R = l -r%L x {X, ® X V X X + X, <g> X X X V ) (21) 

that satisfies ([[]) yields a solution ta of (|20|) ; the general solution is then built 
from and the centraliser of L <g> L. Our theorem has given us the left-hand 
side of (|2lD and a quadratic r-matrix is then given by solving the linear equation 
W c = r ^ F a where pa = ^ + p^)L x /2 and X„X A = F^ x X a . Whereas the 
linear R-matrix involves only Lie algebraic data, the quadratic r-matrix may in- 
volve the group structure through the multiplication X V X\ = F° x X a . Nonethe- 
less, the quadratic r-matrix has been reduced to a linear equation amenable to 
direct solution. 

Finally we mention that for systems obtained by Hamiltonian reduction an 
alternative geometric construction of classical R-matrices exists in terms of 
Dirac brackets. This suggests there is a correspondence between Dirac brackets 
and generalized inverses. This is indeed the case and I will present this elsewhere. 
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